Decay of negative muons bound in ^'Al 
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We present the first measurement of tlie energy spectrum up to 70 MeV of electrons from the 
decay of negative muons after they become bound in '^'^Al atoms. The data were taken with the 
TWIST apparatus at TRIUMF. We find a muon Ufetime of (864.6 ± 1.2) ns, in agreement with 
earlier measurements. The asymmetry of the decay spectrum is consistent with zero, indicating that 
the atomic capture has completely depolarised the muons. The measured momentum spectrum is 
in reasonable agreement with theoretical predictions at the higher energies, but differences around 
the peak of the spectrum indicate the need for 0{a) radiative corrections to the calculations. The 
present measurement is the most precise measurement of the decay spectrum of muons bound to 
any nucleus. 

PACS numbers: 13.35.Bv, 14.60.Ef, 21.30.Fe 



I. INTRODUCTION 

An elementary charged particle can form an atomic 
bound system when it replaces an atomic electron and/or 
a nucleus. Such exotic atoms present interesting systems 
for both basic and applied research, in topics ranging 
from quantum chemistry (e.g. pionic atom chemistry 
[l|]), Coulomb three body systems (positronium ions Qj), 
muon-catalysed nuclear fusion (3]), QED tests (Lamb 
shift [3|), weak interactions (muon capture fsj), strong 
interactions (hadronic shifts 6]), as well as fundamental 
symmetry tests (anti-hydrogen, antiprotonic helium [3]). 

The atomic structure of exotic atoms consisting of a 
heavy negative particle and a nucleus has an unusual 
feature that, because of the larger mass M of the neg- 
ative particle, its characteristic distance scale is smaller 
by me/M than that of the ordinary atoms, and the 
Coulomb interactions are correspondingly larger. The 
average potential energy and the particle velocity, respec- 
tively, are V ~ — (Zq;)^M, and (3 ~ Za, where Z is the 
nuclear charge. Thus, these exotic atoms exhibit bound 
state effects that are significantly more pronounced than 
in ordinary atoms. 

A unique process that takes place in some classes of 
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exotic atoms is disintegration by decay of the short-lived 
constituent. The properties of such short-lived particles, 
such as the lifetime and the decay product energy spec- 
trum, are modified in the presence of the external fields 
of its binding partner. Recently, universal bound state 
principles based on gauge symmetry have attracted in- 
terest. These connect, for example, decay properties of 
electromagnetically-bound exotic atomic states to those 
of heavy mesons, quark-antiquark systems bound by the 
quantum chromodynamic gauge force [^, The 
muonic atom is a system in which the decay proper- 
ties can in principle be calculated very precisely, due 
to the purely leptonic nature of muon decay. Com- 
bined with the strongly enhanced Coulomb interaction 
discussed above, it provides a sensitive testing ground 
for our basic understanding of bound state modifications 
of elementary processes. 

Apart from its own interest as an exotic atomic sys- 
tem, there is currently considerable interest in muonic 
aluminium atoms in the context of searches for muon 
conversion to an electron. Two very ambitious propos- 
als, Mu2e at Fermilab [ll| and COMET at J-PARC [13, 
both propose to use muonic aluminium atoms to look 
for extremely rare lepton flavor violating reactions. Elec- 
trons from muon decay in the bound orbit of a muonic 
atom (DIO) near the end point are expected to give the 
single largest source of background; hence they may limit 
the discovery potential for these projects. 

In this paper, we report a measurement of the elec- 
tron energy spectrum from DIO in the muonic aluminium 
atom using the TWIST (TRIUMF Weak Interaction 
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Symmetry Test) spectrometer. Om^ result for the DIO 
spectrum is the first for aluminium and is the first pre- 
cision measurement covering a considerable part of the 
spectrum in any element, dramatically improving our ex- 
perimental knowledge of DIO. We find that our measured 
spectrum is not adequately described by the existing the- 
oretical calculations, implying the need for inclusion of 
radiative corrections. Our results will confront future 
calculations of DIO, including an upcoming one based on 
the bound-state effective field theory approach [l3|, • 

Historically, the interest in the energy spectrum of the 
decay electron for DIO has been focused on details of the 
nuclear field effects. It was first calculated by Porter and 
Primakoff [T5|. They predicted a Doppler broadening 
of the spectrum by taking into account the momentum 
distribution of the muon in its bound state. In subse- 
quent papers 0,113 , [l^ the calculations were expanded 
to include the modifications to the decay rate and its 
dependence on Z. The spectrum was later recalculated 
with fewer approximations p^ . by taking into account 
higher order corrections to the nuclear potential 12011 and 
nuclear recoil. Furthermore, Herzog and Alder [21| esti- 
mated the effects of bremsstrahlung in the nuclear field 
on the decay electron. Another calculation of the energy 
spectrum and the asymmetry was done by Watanabe et 
al. [22, 23] including tables of numerically calculated val- 
ues for various elements. While these calculations treat 
the muon-nuclear interaction in considerable detail, none 
of them have included the 0{a) radiative corrections 
that arise from the muon-electron interaction. 

Despite the physics potential outlined by the authors 
of the calculations described above, to our knowledge, a 
sufficiently accurate measurement of the energy spectrum 
does not exist. Early measurements [l^, IH, [H, H^] could 
confirm the expected features of Doppler broadening and 
the shift of the spectrum towards lower energies. How- 
ever, these are by no means accurate enough to serve as 
tests for the calculations. Later, data for a very limited 
portion of the spectrum were published for various ele- 
ments in the context of p^—e conversion experiments. Due 
to the nature of those experiments, these results for Cu 

m, s [li,!!^, Ti [HI,!!!,!!, Pb [H,!!! and au [H,!!!] 

focus on the high-energy tail of the distribution and can 
not easily be extrapolated towards lower energies due to 
systematic effects and normalisation problems. However, 
modern [i—e conversion experiments rely in their analysis 
on these calculated spectra (mostly [2l| and [l^) as in- 
put for simulations to estimate the expected background 
in their data. 

The TWIST spectrometer was built for a high- 
precision experiment searching for forms of the charged- 
current weak interaction in the decay of positive muons 
that are not described by the Standard Model. To ob- 
tain a decay spectrum for a free, at-rest muon, /i+ are 
stopped in metal targets and the angular and momen- 
tum spectra of decay positrons are measured. Positive 
(i.e. free) muons must be used for such measurements to 
avoid the influence of the stopping target medium. 



When a negative muon comes to rest in matter, it is 
captured by an atom and replaces an outer shell electron. 
It then cascades almost instantly down to the Is level 
by emitting X rays and Auger electrons. Two processes 
compete for the final fate of the muon: capture by the 
nucleus and DIO. The relative size of these two effects de- 
pends strongly on the properties of the atom in which the 
muon is bound, and this can be determined by observing 
the electrons from the decay. In addition, for the muons 
that do decay before capture, the energy spectrum of the 
resulting decay electrons is modified from that of the de- 
cay of free muons: the energy spectrum is shifted slightly 
towards lower energies as the electron must overcome the 
muon's binding energy (« 0.529 MeV in Al). Also, the 
kinematic endpoint limit is Emax ~ 105 MeV when the 
neutrinos carry away no momentum and the electron re- 
coils against the nucleus. These high- momentum decays 
are very rare, as close to the endpoint the spectrum drops 
sharply oc {Emax - Ef . 



II. EXPERIMENTAL SETUP 

A detailed description of the TWIST detector can 
be found in earlier publications (see [131 and references 
therein). Here, an overview will be given with empha- 
sis on components that are of particular interest for this 
analysis. 

The TWIST spectrometer was located in the M13 
beam fine [Hi at TRIUMF. This channel provided a 
beam of negative muons created by a 500 MeV proton 
beam impinging on a Be target. The beam line was 
adjusted to cloud muons (muons generated in the prox- 
imity of the production target) at a momentum of 29.6 
MeV/c and a rate of w 80 Hz. It had a considerable 
contamination with electrons that was largely eliminated 
at the trigger level. Remaining electrons and a small 
amount of beam pions were identified by scintillators in 
the beam line recording the time-of-flight (TOF) and en- 
ergy deposit of individual beam particles. 

The muons were then transported into the centre of 
the detector and stopped in a thin 71 /im aluminium 
target of 99.999% purity where they decayed. The muon 
range was adjusted using a gas degrader with variable 
density (adjustable ratio of He and CO2), controlled by 
a feedback loop using the measured stopping position. 

The tracks of decay electrons were measured with 
two symmetric stacks of 22 high-precision planar drift- 
chamlaers (DCs) [s^ located upstream and downstream 
of the target (Fig. [T|). In addition, a total of 8 multi-wire 
proportional chambers (PCs) were placed at the very up- 
stream and downstream position of the detector to sup- 
port the event reconstruction by providing timing infor- 
mation. The target was surrounded by another 4 PCs 
(PC5/6 upstream and PC7/8 downstream) to enable the 
measurement of the stopping position of each individ- 
ual muon. The target itself served as the cathode foil of 
the two innermost PCs, thus the gas volumes surround- 
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FIG. 1; (Colour online.) Side view of the TWIST spectrome- 
ter planar chambers and support structures. Muons stopped 
in the target foil, which also served as a chamber cathode. The 
spectrometer is symmetric about the target foil, immersed in 
a uniform 2.0 T solenoidal field aligned with the beam axis. 



ing the target were sensitive. The complete detector was 
contained in a superconducting solenoid magnet, provid- 
ing a highly uniform field of 2 Tesla. 



III. ANALYSIS 



Data Set 



The data set considered for this analysis comprises 
58M triggered events containing 32M beam muons of 
which 5M stop in the target. For those, approximately 
3M decays are observed inside the acceptance of the de- 
tector of which around 1.3M remain after quality and 
kinematic fiducial cuts. 



B. Simulation 

The TWIST detector is reproduced in a Monte Carlo 
(MC) simulation to permit the study of the performance 
of the detector and the reconstruction, and to derive cor- 
rections where necessary. The MC contains a detailed 
description of the detector response and the physics pro- 
cesses affecting the decay tracks inside the detector. This 
simulation is implemented in Geant 3.21 [40| and its cor- 
rectness is verified by direct comparisons between sim- 
ulated and real /i"*" data which is available with high 
statistics. In particular the reconstruction efficiency, res- 
olution, and bias and their phase-space dependence have 
to agree. The complete procedure and results are de- 
scribed in more detail in ^37]- The output of the simula- 
tion is in the same format as the real data and subject to 
the same calibration procedures (where applicable) and 
reconstruction. 



A DC consisted of 80 sense wires at 4 mm pitch sur- 
rounded by Mylar cathode foils of approximately 6 /im 
thickness and separated by 4 mm, filled with dimethyl 
ether at atmospheric pressure. Most DCs were paired 
into modules of two (so-called u and v modules) with the 
central foil shared. A relative rotation of the wire planes 
by 90 degrees allowed for the reconstruction of the posi- 
tion of a hit in the perpendicular plane. The PCs were 
of similar design, but their wire planes were equipped 
with twice as many wires and a "faster" gas (80:20 mix- 
ture of CF4 and isobutane at atmospheric pressure) was 
chosen. In addition to the timing information, the time- 
over-threshold was recorded to provide the amplitude of 
the signal, approximately proportional to the energy de- 
posit. All DC and PC wires, and scintillators were read 
out by time-to-digital converters (LeCroy 1877 TDCs) in 
0.5 ns bins from 6 /is before to 10 /iS after a muon passed 
through the trigger scintillator. The space between the 
chambers was filled with a (97:3) mixture of helium and 
nitrogen, with the relative pressure continuously being 
adjusted to avoid tension and bulging of the chamber 
foils. In total, there was approximately 140 mg/cm^ of 
material from the vacuum of the M13 beam line through 
to the centre of the stopping target. 



C. Event Reconstruction 

To reconstruct an event, the DC hits — signal times 
on individual wires — are first grouped based on timing 
information from the PCs, separately for the upstream 
and downstream halves of the detector. A combinato- 
rial, geometric pattern recognition is performed on the 
hits in such a time window to assign groups of hits to a 
potential track candidate. For each track candidate an 
initial helix fit is performed by using only the spatial in- 
formation of the hits, i.e. the crossing position of a pair 
of hit u and v wires in a module along with the mod- 
ule's z position. This gives the starting values for the 
drift-time fit (DTF). Here, the drift-time information of 
each hit is used together with the space-to-time relation 
tables to position each hit inside the chamber volume. 
The DTF is not a simple fit to a geometrical helix. In or- 
der to obtain better resolution and minimise biases, the 
average energy loss of the particle and hence the dimin- 
ishing radius in the magnetic field is taken into account. 
Also, the deflections caused by multiple scattering are 
included in the fit by allowing kinks i41| at the centre 
of modules, weighted by the width of scattering angles 
expected for the amount of crossed material. The DTF 
results in the particle's time, position and momentum 
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vector at the chamber closest to the target that has con- 
tributed a hit. Therefore the track parameters describe 
the particle where it is first seen, and not at the decay 
vertex. Although this difference only has a small effect 
on the decay time, momentum and the projected angle 
of a track, it is systematic and is taken into account in 
the analysis. 

Event and track selection cuts are then applied to re- 
move background from the spectrum while maintaining 
a minimum of bias in the kinematic parameters of the 
decay. A first set of selection criteria ensures that the 
beam particle is a muon and that it stops in the target. 
The scintillators in the beam line measure the TOF and 
pulse height per trigger, allowing discrimination of muons 
from electrons and pions that are present in the beam. 
The range of the muon is determined by requiring that 
the last hit was recorded in PC6, directly in front of the 
target. In addition, pulse height information of PCS and 
PC6 is used to remove muons that stopped in the gas and 
cause larger energy deposits. Finally, the position of the 
last two PC hits gives an estimate of the muon impact 
point, which should be within 3 cm of the detector axis. 
This ensures that the muon stops in the aluminium tar- 
get, and not in the surrounding support material. Decay 
electrons are then considered starting 500 ns after the 
arrival of the muon. Earlier decays can suffer from re- 
construction problems related to the overlap of the track 
ionisation with that of the incoming muon, which would 
cause an upstream-downstream asymmetry. 

The remaining events are decay candidates and the 
subsequent track selections are aimed at removing poten- 
tial background. Specifically, decay tracks are required 
to have the correct charge sign, and their extrapolated 
coordinates at the target must lie within close proximity 
of the estimated muon impact point. In very rare cases, 
when more than one decay candidate is left at this stage, 
the one closer to the muon position is selected. 

After these selections, kinematic fiducial cuts are ap- 
plied in the E — cos 9 spectrum with E being the total 
energy of the electron and the polar angle with respect 
to the detector axis. These cuts serve to remove regions 
of the spectrum in which the reconstruction is known to 
be less reliable, or where the simulation does not repro- 
duce the data with sufficient accuracy. Specifically, limits 
are imposed on both E and cos0, as well as the trans- 
verse momentum pt (i.e. track radius) and longitudinal 
momentum pi. The exact ranges are adjusted depend- 
ing on specific requirements of statistical and systematic 
uncertainties as well as bias considerations. Kinematic 
fiducial cuts are applied with respect to bin centres, not 
individual tracks. 

For both data and simulation, the track reconstruc- 
tion inefficiency within the fiducial region is of the or- 
der of 10"'^, the energy resolution is between 50 and 150 
keV and the energy bias is less than 10 keV. 



D. Spectrum Unfolding 

Due to the various error sources (biases, resolutions) 
and the limited acceptance and efficiency of an experi- 
ment no measured observable represents the "true" phys- 
ical value. The unfolding procedure tries to solve this 
problem and to find the corresponding true distribu- 
tion from a distribution of the measured observable. 
The main assumption is that the probability distribu- 
tion function in the "true" physical parameters can be 
approximated by a histogram with discrete bins. Then 
the relation between the vector x of the true physical 
parameter and the vector y of the measured observable 
can be described by a matrix Mmig which represents the 
mapping from the true value to the measured one. This 
matrix, usually obtained from a simulation, is called the 
migration (or response) matrix with 

y^Mmig-x. (1) 

In our case the x and y vectors contain the particle energy 
E and polar angle as cos 0. 

The goal of the unfolding procedure is to determine 
a transformation for the measurement to obtain the ex- 
pected values for x using relation ([1]). The most simple 
and obvious solution is the inversion of the matrix. How- 
ever, this method often provides unstable results. Cor- 
relations between bins and statistical fluctuations cause 
the result to be dominated by large variances and strong 
negative correlations between neighbouring bins. 

In the method employed for this analysis [i^l, the un- 
folding is performed by the calculation of the unfolding 
matrix A4^f° in an iterative way that is used instead 
of Ai^^ig- Here is a two-dimensional matrix that 

transforms a vector from the measurement space to the 
space of "true" values. As explained in detail in Ref. [i^ . 
calculating the unfolding matrix this way avoids instabil- 
ities while no assumptions on the migration matrix or the 
shape of the spectrum have to be made. The method has 
been validated by unfolding simulated data. 

To determine A4mig, the MC is used to generate a 
sample of 200M electrons with a decay spectrum flat in 
both COS0 from -1 to 1 and E from to 90 MeV. Since 
the statistical uncertainty of the MC is propagated into 
the flnal result, the sample size is chosen large enough to 
keep this contribution at an acceptable level. 



IV. RESULTS 

A. Lifetime 

The /i~ lifetime is obtained from a maximum likeli- 
hood exponential fit to the track times from the helix 
reconstruction. As noted above, like all other track pa- 
rameters, this time characterises the electron track at 
the first DC encountered. In order to avoid large TOF 
corrections, all tracks that do not have a hit in the first 
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chamber adjacent to the target are removed. As this is 
typically caused by the pattern recognition, it does not 
represent a time-dependent bias. This way only a small 
(sub-ns), angle-dependent correction needs to be applied 
to account for the approximate TOF from the target to 
the first DC. 
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FIG. 2: Reconstructed muon decay time spectrum. The decay 
time binning is 10 ns, the fit range 500 to 7500 ns (error bars 
not shown). 

The measured decay time spectrum and the lifetime fit 
are shown in Fig. [2l The fit range is from 500 ns (below 
which the reconstruction is not guaranteed to be time in- 
dependent) to 7500 ns. In addition to the uncertainties of 
the fit (0.8 ns), two sources of systematic uncertainties 
for this measurement are considered. Detector effects, 
including the non-linearity of the TDCs can be derived 
from a lifetime measurement of free /i"*", and contribute 
with an uncertainty 0.5 ns. Secondly, the muon selec- 
tion uncertainties using a variation on the PC5/6 cuts 
contribute with an uncertainty of 0.7 ns. We then find 

r^i = (864.6 ± 0.8(stat.) ± 0.9(syst.)) ns (2) 

as the lifetime for our muon sample. 

The observed decay time spectrum fits an exponential 
very well, and our result is in excellent agreement with a 
published independent measurement (43j . This indicates 
that the selection of target stops using the PC amplitudes 
is of high purity. 



B. Asymmetry 

The angular asymmetry of the decay spectrum is di- 
rectly proportional to the polarisation of the muon. 



Thus, a measurement of the asymmetry can be used to 
determine the muons' polarisation when the decay oc- 
curs. The TWIST apparatus and analysis have demon- 
strated the capability of measuring the asymmetry of a 
decay spectrum with an accuracy of 1.7 x 10~^ or better 
in the /i+ analysis \^ . 

Contrary to almost fully polarised surface muons 
(muons from pions at rest that decay right at the edge 
of the target), cloud muons have a much lower initial 
polarisation in the opposite direction. The exact treat- 
ment of the mechanisms that lead to a depolarisation 
while the muons travel to the target is not trivial. It 
is assumed that until they come to rest in matter, spin- 
changing interactions (mostly scattering processes) are 
independent of charge [i^. Thus, we can assume that 
the /i^ have a polarisation of about -0.25 before the 
atomic capture. The cascade following the capture will 
depolarise the muons further [i^, \^ . The distribution 
of atomic states in which the capture occurs will depend 
on the atomic structure and is not well calculated. Spin- 
fiip transitions during the cascade produce a large, nearly 
complete depolarisation in the ground state. 

A statistically significant difference in the number of 
decays upstream and downstream is not observed. How- 
ever, a quantitative limit on the residual polarisation can 
be obtained from a fit of the complete angular depen- 
dence of the spectrum and a comparison with the same 



quantity for a data set. With P„ 



'1 and the po- 



larisation depending linearly on the asymmetry, the ratio 
of asymmetries can be used to calculate P^- . Since the 
high-energy region of the spectrum is most sensitive to 
the asymmetry, a range of 31.5 to 52.5 MeV is used here. 
In addition, a small modification for the bound n~ decay 
has to be considered for this comparison; Ref [l^ cal- 
culates that the expected integrated asymmetry should 
be « 7% larger than for a free muon in the considered 
energy range. After this correction we obtain 



P„- 



-0.005 ± 0.003(stat.) ± 0.0017(syst.) (3) 



as the residual polarisation when the muons decay after 
atomic capture. 



C. Energy Spectrum 

Finally, the decay energy spectrum is obtained from 
the unfolding procedure described above. As can be ex- 
pected from the high accuracy and efficiency of the track 
reconstruction of the spectrometer, the unfolding cor- 
rects the spectrum only marginally. Comparing the raw 
and corrected spectra, the most significant difference is a 
small, angular-dependent shift of the energy scale. This 
is to be expected as the measured energy of each track 
at the first DC has to be migrated to its energy at the 
decay vertex. 

The unfolding is performed in the E-cos9 space, but 
since no significant angular dependence of the spectrum is 
observed {P^- ~ at the time of decay), only projections 
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FIG. 3: (Colour online.) Decay energy spectrum. For com- 
parison, the appropriately normalised TWIST spectrum 
is shown as well as the spectrum calculated by Watanabe et 
al. [2^. Most error bars are smaller than the symbols and 
are not displayed. The inset shows the high energy tail of 
the spectrum on a logarithmic scale {with error bars). The 
numerical values are tabulated in Table ID An overall energy 
scale error of 2 x 10~^ is not included. 



are shown, including appropriate acceptance corrections. 
The kinematic fiducial region of 0.54 < | cos < 0.92, 
pi > 14.0 MeV/c, 11.0 MeV/c < pt < 38.0 MeV/c and 
17.5 < E < 73.5 MeV is used. This is slightly different 
from the region chosen for previous analyses (e.g. [37| ) 
in order to optimise the stability of the results towards 
both the low and high energy tails of the spectrum. As a 
systematic uncertainty we assign an overall energy scale 
error of 2 X 10-3(not included in Fig. [3] and Table [H to 
account for an uncorrected difference in the energy scale 
of data and the simulation at the level of a few keV and 
the remaining systematic uncertainties due to misalign- 
ments, differences in resolutions and biases between data 
and MC. 

In Fig. [3] the resulting energy spectrum is shown in 
comparison with the theoretical spectrum calculated by 
Watanabe et al. [2^. The relative normalisation is with 
respect to the number of observed decay events in the 
energy range of 17.5 to 73.5 MeV. In the high-energy 
region (> 60 MeV) the spectrum decreases exponen- 
tially and the observed statistics are very small, requiring 
larger bin widths. To account for this appropriately, the 
abscissa to which a bin is assigned is not the bin centre, 
but given by the integral of the approximate exponential 
probability density function of the population in that bin, 
normalised by the bin width itself (see for example [48|). 



This results in a shift of the energies of the last four bins 
towards slightly smaller values. 
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FIG. 4; (Colour online.) Relative spectrum differences. The 
difference between the measured and calculated spectrum is 
normalised to the theoretical spectrum. The 0[a) radia- 
tive corrections 'is?] are normalised to the free muon decay 
spectrum with an endpoint of ~ 52.8 MeV. 

We see a good agreement in the high energy tail, while 
in the peak region there are significant differences be- 
tween the spectra. At the same time, in the lower end of 
the spectrum our measurement seems to be consistently 
above the prediction. A comparison of these differences 
with the 0{a) radiative corrections (R.C.) for the decay 
of free muons [i^ is shown in Fig. 21 This demonstrates 
that in the region where the differences between the de- 
cay spectra of free and bound muons are small (^ 46 
MeV), the free decay R.C. are applicable, as would be 
expected. To what extent the observed differences above 
the free decay endpoint can be attributed to approxi- 
mations in Watanabe's spectrum, or the lack of radiative 
corrections in the calculation will have to await dedicated 
calculations of such effects for the decay of bound muons. 



V. SUMMARY 

We have measured the decay properties of muons 
bound in ^^Al. The excellent agreement of our mea- 
sured lifetime with published data demonstrates our suc- 
cess in selecting a very pure sample of decays occurring 
within the thin aluminium target material. The depolar- 
isation of the muons during capture must be complete as 
we find negligible angular asymmetry in the decay spec- 
trum. According to our knowledge, the present measure- 
ment of the decay electron energy spectrum is the first 
to be performed on aluminium and the most accurate of 
all measured DIG spectra. Thus this is the first mea- 
surement that can be used to benchmark the calculated 
spectra with precision. In our comparison with the the- 
oretical spectrum (22j we do find differences suggesting 
that 0{a) radiative corrections must be included before 
the assumptions about the basic physics of exotic bound 
systems can be tested. 
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E [MeV] 


events/MeV 


18 


(4.43 ±0.14) X 10^ 


19 


(5.169 ±0.086) X 10^ 


20 


(5.391 ± 0.067) X lO-* 


21 


(5.863 ± 0.054) x 10* 


22 


(6.517 ± 0.060) X 10* 


23 


(6.692 ± 0.054) x 10* 


24 


(7.203 ± 0.051) X 10* 


25 


(7.674 ± 0.053) x 10* 


26 


(7.900 ± 0.053) X 10* 


27 


(8.341 ± 0.050) X 10* 


28 


(8.851 ± 0.056) X 10* 


29 


(9.258 ± 0.058) x 10* 


30 


(9.680 ± 0.055) X 10* 


31 


(1.013 ±0.0061) X 10^ 


32 


(1.052 ± 0.0064) X lO-"' 


33 


(1.077 ±0.0060) X 10^ 


34 


(1.128 ±0.0066) X 10^ 


35 


(1.171 ±0.0068) X 10^ 


36 


(1.193 ±0.0065) X 10^ 


37 


(1.224 ± 0.0070) X 10^ 


38 


(1.262 ±0.0072) X 10^ 


39 


(1.293 ± 0.0069) X 10^ 


40 


(1.301 ± 0.0073) X 10^ 


41 


(1.326 ±0.0075) X 10^ 


42 


(1.371 ± 0.0072) X lO'^ 


43 


(1.381 ±0.0077) X 10^ 


44 


(1.393 ± 0.0078) X 10^ 


45 


(1.401 ±0.0073) X 10^ 


46 


(1.382 ± 0.0079) X 10^ 


47 


(1.352 ± 0.0082) X 10^ 


48 


(1.334 ±0.0082) X 10^ 


49 


(1.264 ± 0.0087) X 10^ 


50 


(1.148 ±0.0084) X lO'^ 


51 


(9.892 ± 0.075) x 10* 


52 


(7.184 ± 0.068) X 10* 


53 


(5.022 ± 0.060) X 10* 


54 


(2.861 ± 0.043) X 10* 


55 


(1.484 ± 0.032) X 10* 


56 


(7.738 ± 0.24) X 10^ 


57 


(4.599 ±0.18) X 10^ 


58.75 (59) 


(1.167 ±0.064) X 10^ 


61.77 (62) 


180 ± 28 


64.79 (65) 


32± 12 


69.39 (70) 


6.7 ±5.0 



TABLE I: TWIST /x"(^''Al) decay-in-orbit energy spectrum. 
The energy coordinates for the last four bins are shifted and 

the bin centres arc given in brackets (see text). An overall 
energy scale error of 2 x 10"^ is not included. 



